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Abstract: A model for the parton distributions in hadrons is derived from sim- 
ple physical arguments leading to an analytical expression for the valence distribu- 
tions. The sea parton distributions arise mainly from pions in hadronic fluctuations. 
Evolving from a low scale with DGLAP gives the proton structure function 
F 2 (x, Q 2 ) in good agreement with DIS data. An extension of the model to describe 
intrinsic charm quarks in the proton is discussed. 

1 Introduction 

When calculating cross sections for hard processes with incoming hadrons, the structure of the 
hadrons is described by parton distributions. The dependence of the parton distributions on 
the hard scale Q of the interaction is described by perturbative QCD evolution. However, their 
dependence on the momentum fraction x at the lower limit for applying perturbative QCD, 
Qo ~ 0.5 — 2 GeV, cannot be calculated perturbatively and must therefore be obtained by other 
means. 

The easiest solution is to fit a parameterization |l|, ||, ||, e.g. of the form fi(x,Q ) = 
NiX ai {l—x) bi (1+Cj y/x+dix) , to experimental data and then use these functions when calculating 
cross sections. The problem with this approach is that the parameters in these functions have 
no direct physical meaning, making it difficult to interpret the results. 

To gain understanding of non-perturbative QCD and the structure of hadrons we have 
developed a physical model [|], §] for the parton distributions at Q . The model is here briefly 
described together with our latest developments. 

2 The model 

The basic assumption in the model is that the probe resolving the hadron has much higher 
resolution than the size of the hadron. The probe will then see free quarks and gluons in 
quantum fluctuations of the hadron. As illustrated in Fig. |I|a, the measuring time is short 
compared to the life-time of the fluctuation which is given by the confinement scale. This 
makes it possible to describe the formation of the fluctuation independently of the measuring 
process. 

2 In proceedings 'Monte Carlo generators for HERA physics', DESY-PROC-1999-02, www.desy.de/~heramc 
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Figure 1: (a) Time scales in probing a parton in a hadron. Kinematics and four- vectors for 
scattering on (a) a valence parton and (b) a sea parton in a hadronic fluctuation. 

Our approach only intends to provide the probability distribution of the four-momentum 
k of a single probed parton. All other information in the hadron wave function is neglected, 
treating the other partons collectively as a remnant with four-momentum r, see Fig. |l]b. The 
parton distribution is described using the light-cone momentum fraction x = ^ which the 
parton has in the initial hadron. Since x is invariant under boosts in the z direction, the same 
will be true for the calculated parton distributions. 

It is most convenient to describe the fluctuations in the hadron rest frame where there is no 
preferred direction and hence spherical symmetry. The momentum distribution of the partons 
should therefore be spherically symmetric. The shape of the momentum distribution should 
be close to a Gaussian as a result of the many interactions binding the parton in the hadron. 
The typical width of this distribution should be a few hundred MeV from the Heisenberg 
uncertainty relation applied to the hadron size. The Gaussian momentum distribution also has 
phenomenological support. The Fermi motion in the proton provides the 'primordial transverse 
momentum', which has been extracted from DIS data and found to be well described by a 
Gaussian distribution of a few hundred MeV width 0. 

2.1 Valence partons 

The probability distribution in momentum space for finding one parton i with mass m; is 
therefore give by the Gaussian 



The Gaussian in the energy component comes from the limited time the parton can be 'free'. 

The parton distributions must fulfill a number of constraints. The normalization for valence 
quarks is given by the sum rules and for the gluons by the momentum sum rule, i.e. 



(fc -m i ) 2 + fc2 + fc 2 +fc 2 



where a 





and 
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There are also the kinematic constraints m 2 < j 2 < W 2 and r 2 > ^Zitn 2 given by the final 
parton j being on-shell or time-like and the remnant having to include the remaining partons. 
This also leads to < x < 1. 

The parton model requires that W is well above the resonance region and that the resolution 
of the probe is much larger than the size of the hadron, i.e. W 3> m p and Qo ^> a^. The 
scale of the probe must also be large enough, Q Aqcd, for perturbative QCD to describe 
the evolution of the parton distributions from the starting scale Qq. 

In OH we integrated eq. [I] numerically to find the parton distributions since the kinematic 
constraints are quite complicated in general. The problem is much simpler if the transverse 
momenta and the masses of the partons are neglected. It is then possible to derive J7| an 
analytical expression for the parton distributions, 

fi(x) = N'(ai) exp ( J erf ( ^-^) where a t = — = « — . (3) 

\ Aaf J V 2^ J m h d h m h m h 

The valence parton distributions for hadrons are here determined simply by the mass and size 
of the hadron! 

The resulting valence distributions for the proton (a w 0.15) and the pion (a ~ 1) are very 
reasonable as shown in Fig. |2|. Note that the pion distributions are very similar to xf(x) = 
2x(l — x) and that one third of the pion momentum is carried by gluons. 




2.2 Sea partons 

The sea partons are described by hadronic fluctuations as illustrated in Fig. |l|c. For example, 
the proton fluctuates into |pvr°) + \mt + ) + ... and the scattering is on a valence parton in 
one of the two hadrons, most importantly the pion. The sea quark momentum distribution 
is obtained as the convolution of the valence parton in the pion, which was already obtained 
from the model, with the momentum distribution of the pion in the hadronic fluctuation. The 
latter is assumed to follow from the same model, i.e. using a spherically symmetric Gaussian 
momentum distribution in the proton rest frame. However, the width a n ~ 50 MeV is smaller, 
related to the longer range of pionic strong interactions or the size of the virtual pion cloud 
around a proton. 
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3 Proton structure function and DIS data 



Using these valence and sea parton x distributions at Qo, next-to-leading order DGLAP evolu- 
tion in the CTEQ computer program was applied to obtain the parton distributions at larger 
Q. The proton structure function F 2 (x, Q 2 ) was then calculated and compared with deep in- 
elastic scattering data as illustrated in Fig. [3|. The model curves were obtained || using the full 
numerical integration including effects from masses, transverse momenta and exact kinematic 
constraints rather than with the simplified analytical expression in eq. (|3|). Furthermore, the 
model parameters are here fitted to the data giving Q = 0.85 GeV, a u = 180 MeV, ad = 150 
MeV, a g = 135 MeV, a T = 52 MeV and 7.7% of the proton momentum carried by sea partons 
(i.e. pion fluctuation normalization). These values are in very good agreement with the defi- 
nite expectations of the model. In particular, the values <7j in eq. ([!]) corresponds to the Fermi 
motion in the proton. (Since a applies in each dimension, one obtains a two-dimensional distri- 
bution with (k±) = 2a 2 in agreement with data on the primordial transverse momentum.) One 
may note that that the u-quark distribution has a ~ 20% larger width, which may correspond 
to a correspondingly smaller available region due to the Pauli exclusion principle. 




Figure 3: Data from NMC |10| (left) and ZEUS (right) on the proton structure function 
F 2 (x, Q 2 ) compared to the model with valence partons only (dashed) and including sea partons 
(full). Extracted from Q where larger x, Q 2 regions are examined. 

The model does remarkably well in view of its simplicity and few parameters (6). Of course, 
conventional parton density parameterizations give better fits, probably mainly due to their 
many more parameters (~ 20). The main advantage with our approach is that it provides a 
model that not only gives a good representation of data, but also provides new insights through 
a physical model for the non-perturbative input distributions. 

The main part of the proton structure is determined by the valence distributions, but the 
sea gives an important contribution at small x. For Q 2 < 1 GeV 2 the DIS formalism and 
the partonic interpretation of F 2 are not straightforwardly applicable. Photoproduction and 
resolved photon processes must here be included in the theoretical description of the data. More 
generally, parton distributions can only be used when the hadron is resolved by a sufficiently 
large scale. Hence, for Q 2 < Qq our parton distribution model can only be applied for processes 
with some other hard scale, e.g. jets in photoproduction, and not to obtain the total photon- 
proton cross section or F 2 extracted from it. 
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4 Intrinsic charm 



Following the model described above we have made a detailed investigation [[12] on how it 
can be extended to describe intrinsic charm quarks in the proton. The main results and their 



comparison with the original intrinsic charm model of Brodsky el al. ]I3[ are briefly discussed 
here. Studies on the possibility to observe intrinsic charm in DIS have been presented earlier 

m. ' 



The idea that our model for parton distributions may be applied to describe intrinsic charm 
is obvious, but there are problems of both conceptual and technical nature. Extending the 
model for generating sea quarks through hadronic fluctuations, illustrated in Fig. [Tic, means 
that one should consider fluctuations \p) — > \A^D ) + .... In old-fashioned perturbation theory 
the particles are on-shell and the fluctuation is suppressed by a factor proportional to 1/AE 2 . 
This provides a description with energy fluctuation in the hadron basis (EFH). 

The large mass of the charm quark and charm hadrons means a large fluctuation AE and 
one could argue that the fluctuation is correspondingly short-lived and that the hadron basis is 
not suitable since hadron states are not be formed. One should then picture the charm quark as 
coming directly from the proton, as in Fig. |l|b, with a remnant containing the charm anti-quark 
(or vice versa) to conserve the charm quantum number. This gives a description with energy 
fluctuation in the parton basis (EFP). 

An alternative description is to impose energy-momentum conservation, which means that 
in the hadron (parton) basis the charm meson (quark) is off-shell. The charm baryon (remnant) 
is external to the hard scattering process and should therefore be on-shell. (The remnant can 
be slightly off-shell if one allows it to pick up some energy in the hadronisation process.) This 
provides two model alternatives with energy conservation in the hadron or parton basis (ECH 
and ECP). 

The kinematic constraints are imposed as previously, but they become harder due to the 
large charm mass. In fact, the ECH model is then not kinematically allowed. The results of the 
other model variations are shown in Fig. [|, with and without the kinematic constraints applied. 
In case of fluctuations in a parton basis there is quark- ant i- quark symmetry, xc(x) = xc(x), 
while for hadronic fluctuations this is not the case since the D meson is softer than the A c 
baryon due to its lower mass. 

As shown in Fig. |]a, the different model variations result in somewhat different charm 
distributions. However, if the kinematic constraints are not applied, then they give essentially 
the same x-distributions as demonstrated in Fig. |]b. They are also then quite similar to the 
x-distribution in the original intrinsic charm model where these kinematic constraints are 
not accounted for. The sensitivity of our results on the treatment of the charm quark mass and 
the constraints demonstrate the importance of taking them properly into account to ensure a 
kinematically allowed final state. 

The results of the models have been studied in more detail Jl2] by examining the four- vectors 
in the process. As expected, the incoming quark k 2 is almost on-shell in the energy fluctuation 
schemes, while in the energy conserving scheme it is off-shell and negative. In ECP, the j 2 
distribution peaks at m 2 and falls very quickly for larger values. In EFP and EFH, however, 
j 2 is larger. This may only be an artifact which can be avoided in the model, but it can also 
signal that perturbative gluon emission must be taken into account. If so, also boson-gluon 
fusion should be included for consistency and then perturbative charm production would enter 
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Figure 4: Momentum-weighted distributions xc(x) for intrinsic charm in the proton obtained 
from variations of our model (EFH and EFP with energy fluctuation in hadron and parton 
basis, respectively, and ECP with energy conservation in parton basis) with (left) and without 
(right) taking into account kinematic constraints. Comparison with the original intrinsic charm 
(IC) model of Brodsky et al. (without constraints). All curves are normalized to unity area. 



such that one would have to define a suitable factorization to disentangle the non-perturbative 
intrinsic charm component. 



Our results here and in |L2] are therefore only a first step towards a realistic model for 
intrinsic charm based on our general idea for deriving parton distributions. 



5 Conclusions and outlook 

We have shown that our new model for parton momentum distributions in hadrons is reasonable 
and can reproduce the measured proton structure function. They can therefore be used as an 
alternative to the conventional parameterizations in practical calculations of hard processes and 
in Monte Carlo generators. The advantage of the physical model over the parameterizations is 
that it provides insights on the non-perturbative dynamics embodied in the model. 

The experimentally observed Jll| flavour asymmetric sea with d > u may be accounted for 
by including all pion fluctuations, but the numerical details remain to be investigated. The 
hadronic fluctuations may also be used to develop a model for the dynamics of leading baryon 
production (c/. Fig. |l|c). 

The application of the model to other hadrons than the proton will be presented elsewhere 
0, but the results cannot be as directly and precisely tested as in deep inelastic scattering 
on the nucleon. However, data on the pion structure are available from Drell-Yan processes in 
pion beam experiments. 
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Applying the model to charm in the proton is interesting, but more complicated due to the 
large charm mass. Although reasonable results are obtained, the x-distribution depends on the 
details in the treatment of the charm quark mass effects and the related kinematic constraints 
imposed to ensure an allowed final state. The proton remnant system plays an important role 
in this respect such that its dynamics may be essential for developing a realistic model. 



References 

[1] H.L. Lai et al, Phys. Rev. D 55, 1280 (1997); 
H.L. Lai et al., |hep-ph/ 9903282 



M. Gliick et al, Eur. Phys. J. C 5, 461 (1998) 
A.D. Martin et al, Eur. Phys. J. C 4, 463 (1998) 
A. Edin, G. Ingelman, Phys. Lett. B 432, 402 (1998) 

A. Edin, G. Ingelman, proceedings of DIS99 to appear in Nucl. Phys. B (Proc.Suppl.) 

G. Ingelman et al., Nucl. Phys. B 206, 239 (1982) 

J.J. Aubert et al, EMC collaboration, Phys. Lett. B 119, 233 (1982) 

M. Arneodo et al., EMC collaboration, Phys. Lett. B 149, 415 (1984) 

A. Edin, in preparation 

M. Arneodo et al., Phys. Lett. B 364, 107 (1995) 

M. Derrick et al., ZEUS collaboration, Z. Phys. C 72, 399 (1996) 

M. Arneodo et al., NMC coll., Phys. Lett. B 364, 107 (1995) 

J.C. Peng et al, FNAL E866 collaboration, Phys. Rev. D 58, 92004 (1998) 

K. Torokoff, 'Charm and strange quark contributions to the proton structure', Master 
thesis, Uppsala University, TSL/ISV-99-0204 

S.J. Brodsky, et al., Phys. Lett. B 93, 451 (1980) 

S.J. Brodsky, C. Peterson, Phys. Rev. D 23, 2745 (1981) 

G. Ingelman, L. Jonsson, M. Nyberg, Phys. Rev. D 47, 4872 (1993) 
G. Ingelman, M. Thunman, Z. Phys. C 73, 505 (1997) 



7 



